Low-energy local density of states of the ID Hubbard model 
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We examine the local density of states (DOS) at low energies numerically and analytically for the 
Hubbard model in one dimension. The eigenstates represent separate spin and charge excitations 
with a remarkably rich structure of the local DOS in space and energy. The results predict signatures 
of strongly correlated excitations in the tunneling probability along finite quantum wires, such as 
carbon nanotubes, atomic chains or semiconductor wires in scanning tunneling spectroscopy (STS) 
experiments. However, the detailed signatures can only be partly explained by standard Luttinger 
liquid theory. In particular, we find that the effective boundary exponent can be negative in finite 
wires, which leads to an increase of the local DOS near the edges in contrast to the established 
behavior in the thermodynamic limit. 
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Interacting one-dimensional quantum wires are well 
studied examples of systems in which the Fermi liq- 
uid paradigm of electron-like quasi-particles is known to 
break down. Luttinger liquid theory predicts that strong 
correlations lead to the remarkable phenomenon of sep- 
arate spin- and charge-density waves as the fundamental 
collective excitations in one dimension at low energies [l[ . 
Experimental confirmation for this picture has long been 
controversial, but now there is some evidence for sepa- 
rately dispersing spin and charge resonances for quantum 
wires on semiconductor hetero-structures 0, Q, quasi 
one-dimensional crystals Q and self-organized atomic 
chains Q as a function of momentum. The density of 
states (DOS) has also been analyzed by scanning tunnel- 
ing spectroscopy (STS) in carbon nanotubes |6|48| and 
self-organized atomic chains This immediately in- 
vites the question, if there are characteristic signatures 
from standing waves of separate spin and charge densi- 
ties, that can potentially be detected in locally resolved 
STS experiments in finite wires. Detailed Luttinger liq- 
uid calculations near boundaries exist, which predict cor- 
responding wave-like modulations in the local DOS and 
a sharp reduction of the DOS near boundaries Tol- 14[. 
However, it is far from clear if these signatures are ro- 
bust in realistic lattice systems, since even minimalistic 
models like the Hubbard chain are not perfect Luttinger 
liquids. The two main reasons for possible discrepan- 
cies are that first of all the assumed degeneracy of spin 
and charge modes in the low energy spectrum can never 
be exact and will be lifted by band curvature and other 
effects. Secondly, it is known that strong logarithmic cor- 
rections will arise from a spin umklapp operator which 
is generically present in ID electron systems with SU(2) 
invariant interactions 15- 17^ . 

The most relevant minimalistic lattice model is the 
Hubbard chain 
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which captures the main aspects of interacting one di- 
mensional electron systems. This model shows signatures 
of spin-charge separation in numerical simulations of the 
momentum resolved DOS [Lsj , which is the central quan- 
tity relevant for photoemission experiments. In this pa- 
per we will now analyze the local DOS as a function of 
position and energy which in turn is relevant for STS 
experiments. 

Recently, several numerical density matrix renormal- 
ization group (DMRG) studies considered the local DOS 
for spinless lattice fermion models in one dimension 0- 
21| , where there is no spin and charge separation. There- 
fore, the complications mentioned above do not arise and 
the agreement with theory is close to perfect in that case. 
As we will show here, the local DOS for a spinful elec- 
tron system away from half-filling is much more complex 
and some key features of the Luttinger liquid theory are 
strongly renormalized. 

The observable of interest is the local DOS to tunnel 
an electron with spin up at a certain energy ui into the 
wire, 



p(x,w) = y j \{u a \ V4, a |0) 5(uj-uj a ) 
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Here, the fermion operator tpt x creates a particle at po- 
sition x on top of the groundstate |0) and the sum runs 
over all states \w a ) with one additional electron. 

In the following we will compare analytical calculations 
for the local DOS from Luttinger liquid theory with sim- 
ulations using the numerical DMRG algorithm 22, 23| . 
For the matrix elements in Eq. (|2|) we use a multi-target 
DMRG with a large number of target states in two dif- 
ferent particle number sectors [lj| and keep track of all 
local fermion creation operators. 

In order to identify the separate spin and charge excita- 
tions, we first focus on the energy spectrum Aw = cj q — loo 
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FIG. 2: (Color online) Total DOS for the excitations of the 
first two levels from Fig. [1] 



FIG. 1: (Color online) Energies Au = ui a — too as a function 
of U for excited states with Nf = iVj, + 1 = 31 particles 
on a lattice of length L = 90. The quantum numbers on 
the right indicate the corresponding mode {m a ,m c }. Dashed 
lines correspond to S — 3 /2 states. For larger U, excitations 
from level 3 (3s, 0c) and level 4 also appear in the low energy 
spectrum. 



for finite wires of length L as a function of interaction U 
as shown in Fig. [T] in units of t = 1. The ground state 
0) is assumed to be a filled Fermi sea with total magne- 
tization S z — 0, so that all particle excited states \u> a ) 
have iV-f = N± + 1 particles and total spin z-component 

S Z = 1/2. 

For U = all excitations are described by simple prod- 
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0) is the low- 



ucts of fcrmion operators 

k n+ i)x where k n = nj^j and \ljq) = cl 
est energy particle excitation. Since the spectrum is ap- 
proximately linear Auj ~ vp(k — fci?), a state involving 
several fermion operators (i.e. a multi-particle excita- 
tion) is nearly degenerate with a single particle excita- 
tion at U = if the sum of excited wavenumbcrs is the 
same, which results in quantized fermion levels shown in 
Fig. [I] For example in level 1, the two states cl JO) and 

cl c{ c^ _ JO) both have approximately the same energy 
lo\ — Wq w vpki. In general there are many multi-particle 
states in each level n, but only one single particle state 

°t,n 



|0) carries all the DOS if E7 = 0. 



The situation changes for finite U. Now all states may 
potentially carry spectral weight in the DOS and the 
near-degeneracy of the levels is lifted. According to the 
Luttinger liquid picture the states are now described by 
integer spin and charge quantum numbers {m Sl m c } with 
energies u) mg t7ric = (m s v s + m c Vc) j yj in terms of the spin 
and charge velocities v s < v c [lTLUj. The spectrum in 
Fig. [T] shows a regular spin and charge spacing with quan- 
tum numbers shown on the right. For each spin/charge 
mode {m s ,m c } there can be several states with approx- 
imately the same energy. The number of states in each 
mode is given by the number of ways it can be created by 



boson creation operators 6^ u , i.e. the product of the in- 
teger partitions of m s and m c . For example, for the mode 
0s 2c in Fig. [TJ we have m c = 2=1 + 1 corresponding 
to the two states created by b\ c and (b\ c ) 2 , respectively, 
which indeed have almost the same energy. On the other 
hand, the near degeneracies of spin modes (e.g. 2s 0c) 
arc significantly split in Fig.[T]by well-known logarithmic 
corrections as will be discussed below 15-17j. 

Regarding the DOS it is now interesting to explore how 
the total spectral weight is distributed among the excited 
states. The compact answer from Luttinger liquid the- 
ory is that the summed up DOS in each mode {m s ,m c } 
should be proportional to a powerlaw of the correspond- 
ing energy, but there are no predictions how the DOS is 
distributed among the states within each mode. For the 
lowest energy modes the situation is still simple, since the 
first level corresponds to one spin and one charge state, 
which should have roughly the same DOS of 1 /2 each for 
small U according to theory. However, the numerical re- 
sults in Fig. [2] already demonstrate obvious deviations 
from this prediction, since the charge mode has a much 
larger DOS for small U . The reason for this discrepancy 
comes from the non-linear band curvature, which lifts 
the degeneracy for finite L even at U — 0, so that the 
interaction has to overcome this energy splitting. Indeed 
for cases where the two states are exactly degenerate at 
U = (e.g. in the thermodynamic limit) the spin and 
charge states have comparable DOS even for infinitesi- 
mally small U. As can be seen in Fig. [2] for L = 90 
the charge state dominates for U < 0.5, which would 
imply that for U < lO 1 ^ the band curvature is dom- 
inant over the interaction effects. The next five states 
from the second fermion level in Fig. [5] show a similar 
crossover behavior with U. Nonetheless, features in the 
local DOS will clearly show the interaction effects even 
in the crossover region as we will see below. 

In the 0s 2c and the 2s 0c modes there are two states 
each as expected. However, only one of the states in 
each mode contributes the overwhelming majority of the 
DOS, while the other state would be practically invisible 
in an STS experiment. The reason that some states in a 
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given mode have zero DOS can sometimes be linked to 
exact symmetries, such as the SU(2) symmetry of generic 
Coulomb interactions. In particular, all particle excita- 
tions \u! a ) in the S z = 1 /2 sector are representatives of 
SU(2) multiplets. The lowest energy state |wo) always be- 
longs to a doublet with S = 1/2. Excitations with charge 
bosons b\ c never change the total spin, but excitations 

with spin bosons b\ may generate higher spin values, 
which can be calculated using the commutation rules of 
the non-abelian SU(2)-Kac- Moody algebra [lj|, since the 
spin bosons correspond to the modes of the SU(2) current 
along the z-direction as discussed in the appendix. For 
example the lowest energy S = 3 / 2 state is given by the 
spin boson excitation -^=[(&j s ) 2 — b\ s ]|^o) with m s = 2 
which is plotted as a dashed line in Fig. [T] Such a to- 
tal spin analysis can be performed for a large number 
of excitations (see also appendix), which is useful since 
states with S > V 2 must have zero DOS due to angu- 
lar momentum addition rules. For example in the fifth 
spin mode 5s Oc, there are 7 states, three of which have 
total spin S = 3 /2 and exactly zero DOS. Interestingly, 
three more states carry only very small spectral weight, 
so that only one state dominates for this mode. This 
indicates that the eigenstates remain the same in terms 
of their bosonic expressions even for U 7^ 0. For charge 
modes there is also exactly one state which carries the 
overwhelming weight in each mode, but the DOS of the 
other states is finite and generally increases with U. 

For a complete analysis it is now useful to turn to the 
local DOS 

Pm„,m e { x ) f° r oach mode {m s ,m c }. Bosoniza- 
tion predicts that the local DOS for a given mode is com- 
posed of a uniform and an oscillating product of spin and 
charge contributions [H| 

Pm a ,mX X ) = M 2 [p",m s ( X )Pc,m e ( X ) 

-cos(2fcpx)p° m# (x)p° mc (s)] . (4) 

In what follows we will focus only on the uniform part 
of the local DOS in Eq. (j4|), since the quickly oscillating 
2fcF-terms cannot be easily resolved in an STS exper- 
iment yet. The amplitudes p„ m ( x ) are slowly varying 
and can be predicted by a simple recursive formula for 
spin and charge (y = c, s) separately [3], 
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where luti x ) = a,^ + cos(2kix). 



(6) 



Here we have defined spin and charge exponents a v = 
(1/K U + K v )/A and b v = (1/K U - K v )/4 in terms of the 
respective Luttinger parameters K v for v = c,s. The 
overall prefactor \c x \ 2 cx (sin -ffi) hc+bs in Eq. (gj) does 
not depend on energy and serves as normalization so that 
Pvm=o = 1- It * s straight-forward to see that the recur- 
sive formula results in powcrlaws for the DOS in the bulk 
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FIG. 3: (Color online) The local DOS of the first few modes 
for L = 92 and U = 1. Points are DMRG data and lines 
are theoretical predictions for K c — 0.9081 and K s = 1.16 
adjusted by a shift as indicated by arrows (see text). Top: 
Local DOS for \ujq). The thick line corresponds to 2|c;c| 2 . 
Lower plots: Uniform part of the local DOS for the first few 
charge and spin modes. 
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14, 2JI. In addition, the for- 



mula predicts slow wavelike modulations in the local DOS 
due to the second term in Eq. ([S]), which also survive in 
the thermodynamic limit near the edge (lp| . 

The local DOS from the DMRG data is shown in Fig. [3 
for the first few modes at i/a-filling. For the lowest ex- 
citation \ujq) the oscillating and uniform parts are the 
same and given by the prefactor |c(cc)| 2 . Already at first 
sight it is surprising to see that the local DOS for all 
modes increases slightly near the boundary, while pre- 
vious works have predicted it to decrease according to 
the boundary exponent [Tol 24 . 25 1. Indeed it must be 
emphasized that the local DOS does not fit the theoreti- 
cal prediction. All curves should in principle be fit free, 
up to one overall normalization, since the local DOS of 
all levels follows from Eqs. (UJ)-©, where the Luttinger 
parameters K C (U) and K s = 1 are known from the ther- 
modynamic Bcthc ansatz. However, in Fig.[3]two impor- 
tant adjustments have been made: First the theoretical 
curves were shifted down for the charge modes and up for 
the spin modes in order to fit the numerical data (indi- 
cated by arrows). This adjustment was already observed 
in the crossover of the total DOS in Fig. [2] due to the 
competition of energy scales (band curvature vs. interac- 
tion) as argued above. Secondly, we find that the spin 
Luttinger liquid parameter must be chosen considerably 
larger than unity K s ks 1.16 for all spin modes in order 
to fit the numerical data corresponding to attractive be- 
havior in the spin modes. This is especially surprising 
since the charge Luttinger parameter from Bethe ansatz 
K c = 0.9081 agrees perfectly with the data without any 
finite size adjustments. There are no other adjustable 
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FIG. 4: (Color online) Left: Luttinger parameters K c (line) 
from Bethe ansatz and K s (points) from the fits to p" m= i(x) 
as a function of U for L = 92. Right: Renormalization of K a 
(points) at U = 1 with the system size L compared to Eq. (0. 



parameters in the fits of Fig. [3j except for one overall 
normalization constant. In return, the results show that 
knowing the local DOS from (numerical) experiments, it 
is possible to extract the effective parameters K s and K c 
from the modulations of only the first spin and charge 
excited states according to Eqs. (g|)-([6|). The oscillating 
parts p° can be analyzed analogously and give the same 
results (not shown). 

In Fig. 0] the behavior of the Luttinger parameters 
from the corresponding fits to the local DOS is shown 
as function of interaction and length at ^-filling. The 
charge parameter from Bethe ansatz K c always agrees 
very well with the data without any additional adjust- 
ments. However, the observed spin parameter K s is con- 
siderably larger than K s = 1. Non-abelian bosonization 
predicts K s = 1 for any SU(2) invariant model, but at the 
same time it is known that a marginal irrelevant opera- 
tor causes corrections to the anomalous dimension which 
only vanish logarithmically slowly with 1/ lnX in the 



thermodynamic limit [151 . Il6j . In abelian bosonization 
such a correction can indeed effectively be modeled by a 
renormalizing spin Luttinger parameter 151-4171 l26l. |27| 



K s — 1 oc g, with g 
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ln(L/L ) (7) 



where Lq is non-universal and depends on the model 
and the quantity of interest. As shown in Fig. U such 
a renormalization description is indeed consistent with 
our data for K s . The parameter K s increases with U 
at a given length L, but decreases slowly as the length 
is increased. The parameter K s appears to be the same 
for all spin modes at a given U and L, i.e. independent 
of energy w. The renormalization of K s is very slow, so 
that exponentially large systems are required to observe 
the thermodynamic limit K s — > 1. The fit parameters 
for K s « 1 + 1.98<7 with In Lq « —6 are outside the range 
what would normally be expected for a spin chain model 
171 ]. Therefore, the particular form of the observed cor- 
rections remains a puzzle. 

Nonetheless, the results of the logarithmic correc- 
tions have interesting consequences. In particular, the 
corrections arc so large, that the boundary exponent 
as = (1/K S + l/K c )/2 — 1 may become negative if 



K s + K c < 2K S K C , i.e. K s - 1 > 1 - K c to lowest or- 
der in the correction, which is indeed the case for small 
U < 2i at L = 92. This results in an increase of the local 
DOS for small energies uj and small distances x near the 
boundary, which is described by a weak powerlaw diver- 
gence p oc x " + c w aB . Such a negative boundary expo- 
nent would also explain the recently observed anomalies 
in the boundary behavior for small interactions in func- 
tional renormalization group studies [28| . 

In conclusion, we have analyzed the local DOS of the 
Hubbard model in the low energy regime. Individual 
states can be classified by separate spin- and charge quan- 
tum numbers. We observe that typically only one eigen- 
state contributes the overwhelming majority of the local 
DOS in each spin/charge mode. The spin and charge 
Luttinger parameters K s and K c can be extracted from 
the modulations in the local DOS of individual excited 
states. While the charge parameter K c agrees well with 
the Bethe ansatz, the spin Luttinger liquid parameter 
is attractive K s > 1 due to large finite size corrections, 
which can only be neglected for exponentially long chain 
lengths. In fact, the corrections to K s are unexpectedly 
strong and may even lead to negative boundary expo- 
nents for moderate interactions U < 2i. The common 
assumption that it is possible to generically use K s = 1 
due to SU(2) invariance is certainly not justified for the 
local DOS. In particular, for finite wires on conducting 
substrates the interactions may be reduced by screening, 
so that the charge Luttinger liquid parameter K c may be 
close to one, while the spin Luttinger liquid parameters 
K s can already be significantly increased, which leads 
to a negative boundary exponent. This would have quite 
dramatic consequences, since the DOS near the boundary 
determines the tunneling between connected wires and 
the renormalization of the conductivity through impuri- 
ties 25[ , which will show an increase at low temperatures 
in this scenario. 
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Appendix: Non-abelian Bosonization 

In order to determine the total spin of an excitation 
it is useful to use non-abelian bosonization in the spin 



channel [15j. In this case the excitations are created by 
the modes of SU(2) currents with a — x,y, z obeying 
the Kac-Moody algebra 



r ra jb 



ie abc J; 



m+n 



(A.8) 



The ground state is characterized by J^^o) =0, Vm < 
0. The total spin operator is given in terms of the m = 
currents 
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(A.9) 



where J ± = J x ± U y . The current modes in the z- 
direction are related to the abelian spin bosons above 

b y J m = \/r b m,s and J-rn = Vt^W for m > °- lt 

is therefore straight-forward to consider the total spin 
of any bosonic spin and charge excitation by using the 
Kac-Moody commutation relations. The charge bosons 
commute with the total spin operator S 2 . Spin ex- 
citations are created by products of spin creation op- 
erators 6^ s acting on |cjo}- The corresponding nor- 
malized states can be labelled by the set of which 
bosons were created {mi, m2, 7773, ....}), e.g. |{3, 1, 1}) = 
-^bl s (b\ s ) 2 \ujo}. Therefore, the matrix elements of 

({mi, 77i 2, m3, ■•■}|S l2 |{ m ii m 2; TO 3i ■•■}) between any two 
such excitations can be evaluated uniquely by the Kac 
Moody algebra (|A.8|) . The J§ operators commute with 
all excitations and the ground state is characterized by 
Jq \ljq) = S z \ujo) = T}\ujo) in our case. For the J ± opera- 
tors we use the Kac-Moody relation in Eq. (|A.8j) with the 
help of computer algebra in order to successively com- 
mute them to the right until the action on the ground 
state is known. This results in a non-diagonal matrix for 
S 2 for each spin mode separately, which can be brought 
into diagonal form. The resulting eigenstates and eigen- 
values are as follows: 
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for the case that the state |wo) has total spin of S = 1/2. 

It is also possible to consider a ojq) state with S — 0. In that case the eigenstates are given by: 
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